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Total correlations of the diagonal ensemble herald the many-body localization transition 
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The intriguing phenomenon of many-body localization (MBL) has attracted significant interest recently, but a 
complete characterization is still lacking. In this work we introduce the total correlations, a concept from quan¬ 
tum information theory capturing multi-partite correlations, to the study of this phenomenon. We demonstrate 
that the total correlations of the diagonal ensemble provides a meaningful diagnostic tool to pin-down, probe, 
and better understand the MBL transition and ergodicity breaking in quantum systems. In particular, we show 
that the total correlations has sub-linear dependence on the system size in delocalized, ergodic phases, whereas 
we find that it scales extensively in the localized phase developing a pronounced peak at the transition. We 
exemplify the power of our approach by means of an exact diagonalization study of a Heisenberg spin chain in 
a disordered field. By a finite size scaling analysis of the peak position and crossover point from log to linear 
scaling we collect evidence that ergodicity is broken before the MBL transition in this model. 


The simple paradigmatic model of a particle hopping on a 
lattice in the presence of disorder significantly advanced our 
understanding of condensed matter systems. It lead to the 
insight that a static disordered potential can lead to a com¬ 
plete absence of diffusion and hence conductance in an iso¬ 
lated quantum system. This is known as Anderson localiza¬ 
tion following its inception by Anderson HI more than half a 
century ago m- The original formulation focused primarily 
on non-interacting systems and in the years following Ander¬ 
son’s work a complete picture was formed: It is now known 
that non-interacting systems in one and two dimensions are 
localized for arbitrary disorder Ei). Anderson also conjec¬ 
tured that a closed system of interacting particles with suf¬ 
ficiently strong disorder would likewise localize and fail to 
equilibrate. This conjecture was only recently put on a firmer 
theoretical footing in a seminal work by Basko, Aleiner and 
Altshuler 0. This has led to a surge in interest in this phe¬ 
nomenon now known as many-body localization (MBL). 

The concept of MBL has been confirmed by a number 
of studies demonstrating that interacting systems 

can display a novel dynamical phase transition between a so 
called ergodic and a many-body localized phase. The MBL 
phase is characterized by robust states protected by the exten¬ 
sively many (approximately) local integrals of motion which 
emerge Many features of this MBL phase have 

since been explored. Lor instance, it has been shown that in 
the MBL phase energy eigenstates typically have low entan¬ 
glement entropy with respect to any bipartition, i.e., satisfy 
what is called an area law ll20l - l^ . This is in stark contrast 
to generic ergodic phases in which the entanglement entropy 
of eigenstates in the bulk of the spectrum exhibits an exten¬ 
sive volume law scaling. Lor an initial pure product states, 
it has also been observed that in many-body localized sys¬ 
tems, bipartite entanglement between two sectors of the sys¬ 
tem grows only logarithmically in time iflOl I23lj29l until an 


extensive value is reached. This differs notably from the usual 
power-law growth in ergodic systems, but also with the non¬ 
interacting case, in which a saturation to a constant is ob¬ 
served. At the same time, many features of MBL are still 
unexplored and their broader connections unknown. 

In this work we go significantly beyond the previous ap¬ 
proaches by applying a powerful and sensitive correlation 
measure to pin down and study the MBL transition. Our fo¬ 
cus is on the time-averaged, dephased states that emerge from 
product initial states once the hopping part of the Hamiltonian 
is abruptly switched on. While fingerprints of the MBL tran¬ 
sition are expected in the correlations of this dephased state, 
their utility depends strongly on the type of correlations con¬ 
sidered. While the behavior of bipartite entanglement is a 
commonly used tool for characterization of phases by the con¬ 
densed matter community, we go beyond this by employing a 
multi-partite correlation measure for mixed states. Quantum 
information theory classifies correlations in quantum states 
as classical correlations, entanglement, quantum correlations 
and total correlations, all of which have distinct physical in¬ 
terpretations and expose subtly different properties ll^ [311 . 
Since we expect the inherently multipartite nature of correla¬ 
tions to play a role in the MBL transition, we argue that the 
total correlations of the dephased state is both a meaningful 
and insightful quantity to investigate it. 

Based on a precise condition for ergodicity we show that 
the total correlations in the dephased state exhibits a differ¬ 
ent scaling with the system size in ergodic and non-ergodic 
phases. In particular, in the disordered Heisenberg spin 1/2 
chain studied we find that the total correlations grow only log¬ 
arithmically with the system size in the ergodic phase, while 
in the MBL phase the growth is linear. Studying via exact 
diagonalization the total correlations averaged over disorder 
realizations and pure product initial states, we show that in 
the crossover region between these two behaviors it develops 
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a pronounced peak with a power-law decay with the disorder 
strength on either side, a key signature to identify and char¬ 
acterize the transition. As a side remark, no peak is visible 
if instead of the total correlations we use the mutual informa¬ 
tion, a measure of bipartite correlations, between the left and 


right half of the system (see also Fig. 10 in the Supplemental 


Material), demonstrating that multi-partite correlations play 
an important role. 

Additionally our study connects the problem of MBL with 
recent research on equilibration in coherently evolving quan¬ 
tum systems. In the past decade this topic has seen an un¬ 
precedented revival of interest mainly due to spectacular ex¬ 
perimental advances in cold atomic physics ||32]|^. In this 
platform, the coherent dynamics can be followed over long 
time scales. In fact, there is evidence for the first experimen¬ 
tal realization of an MBL phase using cold atoms on optical 
lattices m, adding further relevance to the work here. 

Total correlations. The MBL transition has been investi¬ 
gated with a variety of tools, from transport coefficients to 
level statistics. A first diagnostic tool to capture real space 
correlations in quantum states is the growth of the entangle¬ 
ment entropy in the evolution of a product initial state ll24llZ7l . 
Here, in view of the multipartite nature of correlations in in¬ 
teracting many-body systems, we sharpen this approach by 
employing the total correlations T l30l|3l]|35l. In order to 
define T, we first introduce the relative entropy between two 
states p and a defined by S{p\\a) := —tr(p log 2 cr) — S{p), 
where S{p) := —tr(p log 2 p) is the von Neumann entropy. It 
is the quantum analogue of the Kullback-Leibler divergence 
and a very stringent measure of the distinguishability of two 
quantum states ll^ via a result known as quantum Stein’s 
lemma. While not itself a metric, it still upper bounds the 
trace distance via Pinsker’s inequality S'(pllcr) > ||p — cr||f/2 
1361, which captures the optimal distinguishability of quan¬ 
tum states with a single measurement. 

We now introduce the total correlations T: Let V be the set 
of all product states of a A^-partite quantum system, i.e., for 
spin systems, states of the form tt = tti 0 7r2 • • • (8) ttjv (and the 
corresponding analogues for fermionic and bosonic systems). 
The total correlations are then defined as the minimum relative 
entropy between the state and any product state, i.e.. 


T{p) := minS'(p||7r). 
■k£V 


( 1 ) 


It turns out that the unique product state which minimizes the 
relative entropy in the above definition is the product of the re¬ 
duced states Pm obtained from p by tracing out all sites but the 
m-th , i.e., tt = <S)m=iPm EOl- This allows us to compute the 
quantity straightforwardly by making partial traces over the 
partitions of interest. The expression for the total correlations 
becomes 

N 

T{p)='£s{pm)-S{p). (2) 

m—1 

It is useful to point out that for N = 2 the total correlations 
is equal to the mutual information, which has the operational 


interpretation as the work required to erase the correlations in 
p lE). If p is a pure bipartite state, then the mutual infor¬ 
mation is equal to twice the entanglement entropy of p, i.e, 
T(p) = S{pi) + S{p 2 ). We note that although the total cor¬ 
relations defined by Eq. (|^ contains a contribution from the 
diagonal entropy studied in a it also contains contributions 
from all marginal entropies and unlike the diagonal entropy is 
explitly related to the geometric picture of correlations in the 
state under investigation. 

Quantum ergodicity, the diagonal ensemble and many-body 
localization. Leaving aside the problem of a proper defini¬ 
tion of MBL, we take the complementary approach and start 
by defining a property that is a condition for rightfully calling 
a system ergodic. The ergodic hypothesis in classical statis¬ 
tical physics states that ergodic systems explore their phase 
space uniformly such that the infinite time average and the 
microcanonical average should agree (making this precise is a 
subtle issue 1^ ). In quantum mechanics the time and the mi¬ 
crocanonical average can agree exactly only for states that are 
evenly weighted coherent superpositions of all eigenstates in 
a microcanonical subspace lf39l . Hence, we require less and, 
informally speaking, take the standpoint that to call a system, 
i.e., a pair of Hamiltonian and initial state, ergodic (as oppose 
to many body localized) it should explore at least a constant 
fraction of the available Hilbert space. 

Let us now turn this intuition into a clear cut definition. The 
first step is to quantify the explored Hilbert space, we will 
do this based on the dephased or time-averaged state w. For 
a fixed initial state p and non-degenerate Hamiltonian H we 
define 

uj-.= y^\Ef){En\p\Erf{E^\= lim - [ df 
„ r^oo T Jo 

( 3 ) 

where \En) are the eigenvectors of H. This is often referred 
to as the diagonal ensemble, as the off-diagonal elements are 
washed away by the time-average. The dephased or time- 
averaged state is the unique state that maximizes the von Neu¬ 
mann entropy given all constants of motion Ho). If the expec¬ 
tation value of an observable equilibrates on average during 
the time evolution of a system, then the equilibrium expecta¬ 
tion value can be computed from it What is more, un¬ 

der mild additional conditions on the Hamiltonian the follow¬ 
ing is true; If the inverse purity l/tr(w^) of the time averaged 
state, also called effective dimension and participation ratio, 
is high, expectation values of all sufficiently local observables 
equilibrate on average during the time evolution even if they 
were initially out of equilibrium Il^[^l42ll43ll . 

The effective dimension, being a measure for the num¬ 
ber of energy eigenstates that significantly contribute to the 
initial state ll42l l43l . can be interpreted as a measure for 
the explored Hilbert space fraction (as can other moments 
of the energy level occupation distribution, like tr for 
q€Z+ HHl). Instead of demanding a large effective dimen¬ 
sion for ergodicity we only demand the weaker property that 
S{uj) > log(l/tr(a;^)) is large enough to call a system er- 
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godic. 

To identify a reasonable notion of being large enough we 
take inspiration from the theory of random states (although it 
is important to stress that we will not actually base any of the 
later calculations or numerics on Haar random states). For a 
fixed Hamiltonian H and randomly chosen unitarily invariant 
initial states p(0) from the Haar measure on a microcanonical 
subspace of dimension d one can show ll44l Eq. (B 6 )] (com¬ 
pare also Refs. Il45lj48l ) that for some C > 0 

Pr(S'(a;) < log2(d/2)) < 4 exp(-C'(i/log 2 (d)^). (4) 

That is, random states typically explore at least half of the 
available Hilbert space in the sense that typically S{uj) > 
log2(d/2). 

For our condition for ergodicity we relax this fraction of 
1/2 to a constant fraction of the available Hilbert space. To 
make this meaningful we have to speak about families of sys¬ 
tems of increasing system size N, specify what we mean by 
available Hilbert space and describe the class of initial states. 
As is common in localization studies we take the subspace 
of dimension d corresponding to a fixed filling p G [ 0 , 1 ] 
or magnetization 2 p — 1 as the available Hilbert space. We 
then consider initial states that are pure product states with 
definite local particle number or magnetization from that sub¬ 
space, which can be thought of as ground states of appropriate 
“easy” Hamiltonians. We say that a family of such systems 
should be considered ergodic only if most such product ini¬ 
tial states explore at least a constant fraction of the fixed fill¬ 
ing/magnetization subspace in the sense that for some A > 0 it 
holds that S(tv) > log 2 (Ad). Note that this is less restrictive 
than demanding that l/tr(a;) > Ad, as S(tv) > — logtr(w). 
For families of disordered systems we demand that the same 
condition is fulfilled with high probability also with respect to 
the disorder average. 

Scaling of the total correlations. We now turn to demon¬ 
strating that the total correlations in the dephased state can be 
used to pin down and better understand the transition point 
from an ergodic to the MBL phase. The key signature we ex¬ 
ploit is the scaling of T(u;) with the system size N. Inspecting 
Eq. Q one might expect that the total correlations in the de¬ 
phased state T(w) should generally scale extensively with TV, 
i.e, for large N one should have to leading order 

T{uj) oc N, (5) 


For a quantum spin chain of local dimension 2 at half fill¬ 
ing rj = 1/2 the available Hilbert space dimension is d = 
(^ 2 ) = > v^e -2 2 ^/y]V and < 

log 2 2 = 1, SO that one finds at most the logarithmic scaling 

T{uj) < log2(TV)/2 - log 2 (A v^e"^). (7) 

This is what we observe in the ergodic phase of the model we 
consider. One furthermore retains a logarithmic scaling for 
ergodic spin 1/2 systems for all other constant fillings 77 G 
[ 0 , 1 ] if another mild condition is satisfied that can also be 
motivated from ergodicity, namely that for some A' > 0 

N 

S{ojjn) < N s{r]) + X'log2N ( 8 ) 

m—1 


where s{x) := —x log 2 (a;) — (1 — x) log 2 (l — x) is the bi¬ 
nary entropy function. Eq. ([^ says that the sum of the local 
entropies of the time averaged state should not grow much 
faster than one would expect for the given filling fraction rj. 
The generalized Stirling formula implies that 

log 2 d = log 2 >Ns{t])- c{rf), (9) 

with c(?7) = 21 og 2 (e) — log 2 (p) — log 2 (l — rfj. Inserting 
this and Eq. ([^ into (|^ yields a logarithmic scaling with the 
system size: 


T(w) < A'log 2 TV-log 2 (A)-f 0 ( 77 ) (10) 


This sub-extensive scaling can also be understood intu¬ 
itively: The transport present in a ergodic systems correlates 
the different parts of the system to the extent that they ap¬ 
pear, for most times during the evolution, so mixed that the 
distinguishability from the closest product state only grows 
logarithmically. 

Model used for numerics. A model which is known to ex¬ 
hibit a crossover between an ergodic and a MBL phase is the 
Heisenberg spin chain with random field in the z direction [7] . 
The Hamiltonian of this model is given by 


H = 


N 


2=1 


( 11 ) 


as T{uj) involves the sum J2m=i of the N subsystem 

entropies. Indeed, this is the behavior we find in the MBL 
phase of the model we consider below (see also Fig. in the 
Supplemental Material). 

If a family of disordered systems is ergodic however, then 
for some constant A > 0, for most product initial states, and 
with high probability over the disorder average 

N 

T{oj) < ^ S{u}m) - log 2 (Ad). ( 6 ) 

m—1 


where the hi represent identically distributed static fields on 
each site i uniformly distributed in the interval [—h,h]. In 
what follows we adopt periodic boundary conditions and set 
Jz = 1 , so that a family of systems is completely charac¬ 
terized by the XX type coupling constant J and the disor¬ 
der strength h. For all values of the parameters, the model 
conserves the total magnetization Sz along the z direction, so 
in the numerics we have chosen the subspace with Sz = 0 , 
also referred to as half filling, i.e, 77 = 1/2. We take as our 
initial states all j^roduct eigenstates of the on-site part of the 
Hamiltonian from this subspace. We then compute 
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FIG. 1. For high h, T decays as a power-law. The dashed lines 
are fits to the data points with h > 10 yielding exponents of —0.9(1) 
(consistent with the expectation of h~^ corrections). The inset shows 
the position of the peaks from Fig.(see also Fig.|^in the Supple¬ 
mental Material). The extrapolated position of the peak indicates the 
onset of MBL at around h = 3.8. The fact that the different T/N 
curves overlap for large enough h indicates that T scales essentially 
linearly with N in this regime. 



h 

FIG. 2. Subtracting the values of T for h = 0.2 (which is just outside 
the integrable region around h = 0) from T one can see that for 
increasing system sizes the individual curves fall on top of each other 
in an increasingly large region of h values (shaded regions) during 
the approach to the peak. This increase is well captured by a power- 
law with exponent 2.7(2) (dotted line guide to the eye oc 

for each initial state the diagonal ensemble w and T{uj). Av¬ 
eraging over all such initial states and disorder realizations 
yields T{oj). The numerics were performed using standard 
libraries for matrix diagonalization. We use 10, 000 disorder 
realizations for each disorder amplitude \h\ and system size 
N, except for the case of = 16 where 1000 realizations per 
point were computed. 

Results, discussions and conclusions. We concentrate on 
the case J = = 1. For that case the MBL transition in 

the model o was predicted to be he G [2,4] by Huse and 
Pal 121, with the best estimate based on energy resolved calcu¬ 
lations being he = 3.72(6) ifT^ from an analysis of spectral 
statistics. Turning to the results we obtained, in all calcula- 



FIG. 3. The difference of the average total correlations T and the 
best possible linear fit for different values of h. The crossover from a 
nearly linear scaling to a sub-linear scaling clearly happens between 
h = 2.8 and h = 2.4. This result is robust against omitting data 
points for large or small values of N and equally holds for affine 
fits instead of linear ones. This gives us high confidence in this re¬ 
sult, which indicates the existence of an extended non-ergodic but 
not yet many-body localized region.The inset shows the data before 
subtracting and the fits. 


tions performed we observe that the total correlations when 
plotted versus h show an initial growth at low h towards a 
maximum and then decrease monotonically at higher disorder 
with a power-law decay with an exponent of roughly —0.9(2), 
i.e., T (X (see Fig . [^. This is consistent with the 

expectation of h~^ corrections from perturbation theory and 
the behavior in the non-interacting case = 0 (see the Ap¬ 
pendix). The position of the maximum is size dependent, and 
can be extrapolated to be he « 3.8 in the thermodynamic 
limit. This is in excellent agreement with the best known ap¬ 
proximation of the transition available in the literature ifT^ . In 
turn, by rescaling T/N it appears that all curves collapse onto 
a single master curve for h > he (see Fig.j^in the Appendix). 
Since on qualitative grounds we expect the many-body local¬ 
ized phase to be characterized by (i) linear scaling T x N, 
(ii) shrinking localization length as h increases, it is natural 
to identify the MBL transition with the peak (see the inset of 
Fig. B- The analysis for L < he is more complex: scaling and 
data collapse for T / log(W), expected in an ergodic phase, are 
observed only for h < 2.6(2) (see Fig.[^and[^. For low dis¬ 
order (up until h K, 2.0 for the system sizes we can access), 
we see a power-law increase of T with h (see Fig. [^ with an 
exponent of about 2.7(2). While for h < 2 the system is defi¬ 
nitely in an ergodic phase (compare Fig.|^in the Supplemen¬ 
tal Information), the analysis resented in Fig. suggest that 
ergodicity is is broken only around h = 2.6(2). This is con¬ 
sistent with an intermediate extended yet non-ergodic phase 
MM before full MBL sets in. Due to the small system size 
finite size effects cannot be completely ruled out. Neverthe¬ 
less, recent work has demonstrated an intermediate level stat¬ 
ics implying non-ergodic extended states in precisely the 
crossover region indicated by our numerics. We believe that 
our work constitutes evidence of an intermediate non-ergodic 
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region before the onset of MBL. 

Conclusions. The numerical simulations performed to¬ 
gether with our analytical arguments show that the total corre¬ 
lations in the diagonal ensemble signal both ergodicity break¬ 
ing and the MBL transition in a quite spectacular way. In stan¬ 
dard critical systems it is known that the multi-partite correla¬ 
tions of the system rearrange as the system is pushed across an 
equilibrium phase transition ll50l . Undoubtedly the transition 
from an ergodic to a MBL phase is a highly non-equilibrium 
phenomenon which is poorly understood at present. Our ap¬ 
proach exposes how this transition goes along with a reorga¬ 
nization of correlations in the dephased state via significant 
change in scaling with N. We expect this behavior to be 
generic and believe that the methodology outlined here is very 
promising to study MBL and ergodicity breaking phenomena 
in a variety of many-body quantum systems. In a follow up 
study we will investigate the possible multi fractal nature of 
ergodicity breaking in a way inspired by Ref. BTl l49l . 
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APPENDIX 

Total correlations in Anderson localization 

In this section, we complement the results presented in the 
main text by discussing the total correlations for the time av¬ 
eraged state in the random XX model, with Hamiltonian 

N 

+ ( 12 ) 

i=l 

The above model is equivalent to the model in the text for 
Jj = 0, and serves as a simple testbed in which very large 
system sizes can easily be probed. Specifically, for suitable 
boundary conditions, the above model is equivalent to the free 
fermionic model 

H = pMf, (13) 

which we take as the basis for our analysis of non-interacting 
disordered models. Here, / = (/i,..., /at)^ is the collec¬ 
tion of free fermionic annihilation operators of N fermionic 
modes. In the above quadratic form, the kernel M is given by 

M := —2 J circ(0,1,0,..., 0,1) — 2 diag(/ii,..., /iat), 

(14) 

again with {hi,..., drawn uniformly random from 

[—h,h]^. This real symmetric matrix can be diagonalized 
as 

M = ODO'^, O G 0{N), (15) 

with D being real and diagonal. More generally, unitary trans¬ 
formations U G U{N) from one set of fermionic operators to 
another one can be allowed for, and in all what follows, or¬ 
thogonal transformations can be replaced by unitaries. For 
the present purposes, this is unnecessary, however. 

Gaussian states p of systems of massive fermions (no Ma- 
jorana fermions are considered) can be captured in terms of 
correlation matrices C{p) > 0, with entries 

Cip)j,k ■=fT^{f]fkP), j,k G N}. (16) 

For a given correlation matrix C{p) of such a state p the cor¬ 
relation matrix of the state a expressed in the basis in which 
the Hamiltonian is diagonal is given by 

C{a) = 0^ C{p)0. (17) 

It is easy to see that the correlation matrix of the infinite time 
average w of the initial state p represented by C{p) is then 

C{oj) = 01f{0'^ C{p)0)0'^, (18) 

where H is the map that projects a matrix onto its main diag¬ 
onal. 

The ground state correlation matrix can again be expressed 
in terms of the kernel M of the Hamiltonian form, as long as 



FIG. 4. Data points on the right are the averaged total correlations 
of the free model with Jz = 0. The behavior is again well captured 
by a power law with exponent —0.9(1). For comparison on the left 
is the data shown in Fig. [T] of the main text. Dotted lines are power 
law fits. 


{0} ^ spec(M): Then the ground state p is unique and has 
the correlation matrix 

C'(p) = O sgn(diag(M)) (19) 

as again can be verified by expressing the Hamiltonian in the 
appropriate basis. 

Entropies of Gaussian states can be computed from their 
correlation matrices. Making again use of the binary entropy 
function one finds that any such Gaussian state p of N modes 
with correlation matrix 1 > C{p) > 0 has the von-Neumann 
entropy 

S{p) = tr(s(C'(p))), (20) 

as can be seen by exploiting suitable orthogonal mode trans¬ 
formations and the unitary invariance of the von-Neumann en¬ 
tropy on the level of quantum states. Hence, the total correla¬ 
tions of the time averaged state oj are found to be 

T(p) = tr(s(n(0n(0^C(p)0)0^)) 
-tr(s(n(O^C(p)0))). 

That is to say, both entropies of the infinite time averaged state 
and its reductions can be conveniently computed. In this way, 
once the correlation matrix has been identified, the total corre¬ 
lations measure can be immediately obtained. Drawing i.i.d. 
random vectors {hi,..., hj^)'^ uniformly from [—h, h]^ as 
in the main text, one can very clearly identify the power law 
decay of the total correlations. Numerically, system sizes of 
N = 100, can easily be accommodated in this way, finding 
that logT/N log(/i) can be well fitted with an affine func¬ 
tion, reflecting a power law, again with exponent —0.9(1) (see 
Fig.j^. 
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Further plots 

In the following we present some additional plots for the 
Heisenberg spin chain with random held in Eq. o of the 
main text with with J = Jz = 1. 




N 

FIG. 7. Scaling of the averaged total correlations T for the transla¬ 
tion invariant case without disorder h = Q and in the MBL phase at 
h = 4. The dashed lines are linear fit to the /i = 4 data and the graph 
of logj N (no fit) for comparison. As in the ergodic phase we find a 
logarithmic scaling even at h = 0. 


h 


FIG. 5. Semi logarithmic plot of the intensive averaged total corre¬ 
lations for the Hamiltonian in Eq. 0 with J = J 2 = 1 versus h for 
system sizes between 6 and 16. The merging of the different T /N 
curves for large h demonstrates that we can reliably determine that 
T exhibits the linear scaling predicted in Eq. l|^ with the system size 
N down to h ~ 2.6. 



h 


EIG. 6. Semi logarithmic plot of 2^'^ IN versus h for system sizes 
between 6 and 16. The merging of the different curves demonstrates 
that we can reliably determine that T exhibits the logarithmic scaling 
predicted in Eq. 0 in the main text with the system size N up to at 
least h ~ 2.Q. The more elaborate analysis performed in Fig.of 
the main text even pushes this to h — 2.6(2). 



h 

FIG. 8. The standard deviation a of T with respect to both the 
disorder average and the average over product initial states divided 
by T as a function of h. Well within a phases either (MBL or ergodic) 
T is self averaging. However, if closer to the transition the quantity 
is affected by rare events. We expect the region around the peaks, in 
which a/T scales linearly with N, shifts further to the right, into the 
region of the phase transition with increasing system size. 
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l/N 

FIG. 9. The figure shows the position of the peaks in Fig.|^of the 
main text and is an enlarged version of the inset of Fig. [^of the main 
text. 



h 

FIG. 10. Compared to the total correlations, the mutual information 
Mip) ■- S'(p[i,jv/ 2 ]) + S{piM/ 2 +i,N]) - S{p), when computed 
in the time averaged state and averaged in the same way, appears 
to be mostly featureless in the parameter range of the localization 
transition. This supports our point that the multipartite nature of the 
total correlations is the reason for why it signals the transition so 
well. The plot shows an average over 100 disorder realizations. 












